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We show that a spacetime satisfying the linearized vacuum Einstein equations around a type D 
background is generically of type I, and that the splittings of the Principal Null Directions (PNDs) 
and of the degenerate eigenvalue of the Weyl tensor are non analytic functions of the perturbation 
parameter of the metric. This provides a gauge invariant characterization of the effect of the 
perturbation on the underlying geometry, without appealing to differential curvature invariants. 
This is of particular interest for the Schwarzschild solution, for which there are no signatures of 
the even perturbations on the algebraic curvature invariants. We also show that, unlike the general 
case, the unstable even modes of the Schwarzschild naked singularity deforms the Weyl tensor into 
a type II one. 

PACS numbers: 97.60.Lf, 04.20.-q 


I. INTRODUCTION 


Understanding how the background geometry is affected by gravitational waves, here meaning solutions hat to the 
linearized vacuum Einstein’s equation (LEE) around a vacuum solution gab- 

1 . ^ablOcd “1“ e/lcci] 

lim- = 0, 

e^O e 


or, more explicitly, 


V^Vchab + VaVb{g^‘^h,d) - 2V=V(,hb), = 0, 


( 1 ) 


is a rather intricate problem, partly due to the gauge issue of linearized gravity, that is, the fact that given a solution 
hab of (1) and an arbitrary vector field h'^f, defined by 

^ab ~ ^ab “b ^VQab: (2) 


is also solution, although physically equivalent to hab- We are only interested in the equivalence classes of solutions of 
(1) under the equivalence relation ^ hab + £v9ab, and on functionals of hab that are gauge invariant, i.e., depend 
only on the equivalence class of hab- The linear perturbation ST of a tensor field T that is a functional of the metric, 
transforms as ST —>■ ST + TyT under the gauge transformation (2), therefore only constant scalar fields and tensor 
products of the identity map S^ are gauge invariants. 

A possibility explored in [1] is to parametrize the equivalent classes of metric perturbations in terms of the pertur¬ 
bations of a set of curvature scalar fields that vanish in the background. In vacuum, there are only four functionally 
independent algebraic invariants of the Weyl tensor Cabcd and its dual (recall that left and right Hodge duals 

of Cabcd agree), these are 


^_1 ^ ^abcd 

•— 4g^a&C(i^ 

C . 1 cd/^ ef/o ab 

+ •— gQ^ab ^cd ^ef 7 

^ _ 1 ^abcd 

V- ■— J^'^abcd^ 5 

C . 1 cd/^ 

“ 96*^06 '^r.d '-^cj 


( 3 ) 


For the Schwarzschild black hole only Q_ and C- vanish in the background, and its is shown in [1] that SQ- 
parametrizes the space of odd (also called vector, see, e.g., [2] [3]) metric perturbations, and that SC- oc SQ- and 
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therefore adds no information. On the other hand, under even (or scalar, [2] [3]) perturbations, every gauge invariant 
combination of the perturbed scalars in (3) vanishes identically. For this reason reason, the gauge invariant combina¬ 
tion (9M — 4r)dQ+ + 3r^SX, which involves the differential invariant X = (l/720)(VaC'{,cde)(V“C'{,cde) was added to 
5Q- in [1] to parametrize the entire set of perturbations using geometrically meaningful quantities. A natural question 
to ask is what -if any- are the effects of the even perturbations on the curvature itself; in other words, do we really need 
to look at differential invariants to find a geometric signature of the perturbation? That contractions of products of 
the curvature tensor hide vital information in Lorentzian geometry is not a surprise: pp-waves are an extreme example 
of non flat vacuum metrics for which every scalar made out of Caked and arbitrary covariant derivatives of it vanish [4]. 


The Weyl tensor of a generic metric is of type-I in the Petrov classification. This means that the eigenvalue problem 




(4) 


admits three different solutions, with Ai -I- A 2 -I- A 3 = 0 or, equivalently, that the equation 

kieCa]bcidkf]k’’k‘' = 0, k‘^ka = 0, (5) 

admits four solutions spanning four different null lines (called principal null directions or PNDs). Type D spacetimes, 
instead, are characterized by the fact that the eigenvalue equation (4) admits three linearly independent solutions 
with Ai = A 2 = — 5 A 3 , a condition that turns out to be equivalent to the existence of two so called double PNDs, that 
is, two non proportional null vectors satisfying an equation stronger than (5): 

Cabc[dke]k^k^ = 0 , = 0 . ( 6 ) 

It should be stressed that equations (5) and ( 6 ), being homogeneous, do not define (null) tangent vectors at a point 
p of the spacetime (i.e., elements of TpM), but instead one dimensional subspaces through the origin of TpM (this 
is why se talk of null directions), that is, points in the projective space P(TpM). We may think that type I (D) 
spacetimes have four (two) smooth valued functions that assign to every p G M a point in P{TpM). Now suppose 
that gab{e) (e in an open interval around zero) is a mono-parametric family of vacuum solutions of type I for e ^ 0 
and type D for e = 0. The tensor 


hab = lim 
£^0 


gabje) -gafc(O) 
€ 


(7) 


satisfies (1) and the four PNDs of Cabcd{e) coalesce pairwise into two PNDs as e —>■ 0. We can choose solutions 
fc“(e), ^2 (e); ^3 (e) and k%{e) of k\^eCa]bc[d{^)kf]k^k^^ = 0 such that fci(O) = A: 2 ( 0 ) and ^ 3 ( 0 ) = ^ 4 ( 0 ) are the two non¬ 
proportional solutions of the type-D equation Cabc[d{^)ke\k?’k'^ = 0. Similarly we can label the eigenvalues Ai(e), A 2 (e) 
and A 3 (e) of (5) such that Ai(0) = A2(0) = — 4A3(0). This suggests using either the splitting of type D PNDs in 
P{TpM) into a pair of PNDs, or the eigenvalue splitting 


±(A2(e) - Ai(e)) 


( 8 ) 


to measure the distortion of the unperturbed e = 0 metric. Note that the sign ambiguity in ( 8 ) comes from the 1 GG 2 
freedom in labeling the two eigen-bivectors (4) that degenerate in the e = 0 limit. We have found that for e ~ 0, 
these are appropriate, gauge invariant estimates of the distortion of the geometry by a gravitational wave involving 
only the Weyl tensor. 


The reason why we expect these quantities to be non trivial is that the bulk of vacuum solutions of the Einstein’s 
equation are type-I, while algebraic special solutions comprise a zero measure subset A. The curve gabip) will gener- 
ically be transverse to A at gab{0) (as assumed above), its “tangent vector” hab will stick out of A (more on this in 
Section IV C). 

An speciality index S was introduced in [7] based on the fact that a vacuum spacetime is algebrically special if and 
only if = Q^; it was defined as S := C\IQ\ and the departure of this index from unity is regarded as a measure 
of the “degree of non-speciality” of a metric. Although useful in numerical evolutions of the full Einstein equation, 
5 — I vanishes identically in linear theory, suggesting the -wrong- conclusion that linearly perturbed algebraically 
special spacetimes remain special. This error was nicely clarified in the paper [9], taking advantage of the fact that 
a mono-parametric set of vacuum solutions within the Kasner family, with gab{0) being type D, can be explicitly 
constructed. From the results in [9] is to be expected that the PND splitting will be a power series of This 
non analyticity on e comes from the fact that the associated algebraic problem involve radicals that vanish for e = 0 [9]. 
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Stationary electro-vacuum black holes have the algebraic symmetries of type D. Deviations from these metrics 
represent a number of very important astrophysical processes and, as argued above, will not to be algebraically 
special. Thus, in the context of black hole perturbation theory it is important to study whether or not the perturbative 
techniques capture these algebraic aspects of the geometry. In numerical application this question has been analyzed 
in [7],[8]. On the analytical side in [9] for the Kasner type D solutions. In the present work we focus on the 
problem of finding expressions for the eigenvalue and PND splittings of the type D background under arbitrary linear 
perturbations. 

As an application, we consider the effects of linear perturbations of the Schwarzschild black hole solution and show 
that the even or scalar gravitational waves (those that do not leave a visible trace in the perturbed algebraic curvature 
invariants) do affect the Weyl tensor by turning it into type I, and that this effect is accounted for by the proposed 
gauge invariants above and we also show that these invariants mix harmonic modes. As far as we are aware, these 
signatures of the black hole linear perturbation are analyzed in this form for the first time here. 

Kerr black holes and Chandrasekhar algebraic special modes in Schwarzschild naked singularities are also briefly 
considered. 


II. PETROV TYPES 

Although the simplest approach to the Petrov classification is accomplished using spinor methods, perturbation 
theory is much more tractable in tensorial language, which is the one we adopt in this paper. Newman-Penrose 
equations will also be avoided, although complex methods and, in particular, a complex null tetrad will be used. 
For the sake of completeness, and to clarify some aspects of the linearized theory calculations, we briefly review the 
eigenvalue and PNDs approaches to the Petrov classification in this section. 


A. Eigenvalue approach 

The eigenvalue approach to the Petrov classification regards the Weyl tensor as a linear map —>■ C°‘^cdX^‘^ in 

the space of rank two antisymmetric tensors (also called “bivectors”), and analyzes the eigenvalue equation (4). The 
six dimensional bivector space is real isomorphic to the complex three dimensional space of self-dual bivectors (SDB), 
which are those satisfying *Sab '■= = —iS°'^. The isomorphism is given by 

j^ab ^ ^^ab^^j^cd j^ab^ 

its inverse being 

i -h c.c.) (10) 

The space of SDBs, in turn, is isomorphic to u^, the complexification of the space of vectors orthogonal to a given 
unit time-like vector u°'. The isomorphism is 


j^ab 

and its inverse is 

1“'' = -b i€^\au^X^. (12) 

If A°'b is a Lorentz transformation of unit determinant 

= g-\ A%A%A%A‘^,e^’«”* = tabcd, (13) 

one obtains from (13) that (A“^)“t, = Af,“ and 

A%A\eP<^cd = e“^,A”,A^d. (14) 

Equation (14) implies that the unit determinant Lorentz transformations commute with the map (9), as expected 
from 


1 \^ab V _ vab v _ a v'a \r 

-^ab — -A -^ab — '^-A A.q,. 


( 15 ) 
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The second equality above, obtained from (12) together with UaX°' = 0, indicates that A acts on as an S'0(3,C) 
transformation —>■ A^'i^X^. With the help of (12) and (11) we find that this gives an isomorphism between 

SO{3, 1)^ (the group of unit determinant Lorentz transformation preserving time orientation, i.e., A°'t,UaU^ < 0 if ■u'^ 
is time-like) and 5'0(3, C): 


A\ = A% - A“o AOfc + ie°%i A% A'^, (16) 

where a zero down (up) index means contraction with {—ut). To gain some intuition on the isomorphisms (9) and 
(11), note that if Xab is the electromagnetic tensor, then X°' = E‘^ + iB°', the electric and magnetic fields measured 
by an observer moving with velocity If A“f,t 6 ^ = u“, then A is a rotation, and (16) reduces to 


A°' k — 11°' Uk A° k 


(17) 


which belongs to S'0(3,]R) and simply rotates E and B independently. Otherwise, A is a boost and (16) a complex 
rotation mixing E and B, as expected. 

We can replicate the above constructions for the self-dual piece of the Weyl tensor Cabcd by regarding it as an element 
of the symmetric tensor product of bivector space, and using the fact that that left self-duality implies right self 
duality. Define 


Cabcd ■= 5 (Cabcd + Cklcd) , 


and introduce the map Q : ^ 


Q°a := -C°bcd 


(18) 


(19) 


The above equation can again be inverted [5], 

^Cabcd ^^[aQb][d^c] T 9a[cQd]b 9b[cQd]a T 'i^^abefH' d]^ T i^cdefH' ^[aQb]^ ^ 

and the eigenvalue problem (4) is easily seen to be equivalent to jC°^cdX°'^ = \X°^, and also to 


Q°bX^ = \X°, X° = X°°Ub. 


a \ra _ \^ab^ 


( 20 ) 


( 21 ) 


We should stress here that, although Q°b and X^ were defined using a particular unit timelike vector u°, the eigenvalue 
equation (21), being equivalent to (4), gives covariant information and therefore is fully meaningful. Let {e° = 
u°, e°, 62 , 63 } be an orthonormal tetrad {gabC^e^p = diag(— 1 , 1 , 1 , 1 )) adapted to u°, 


a _ '=^0 


k° = 


a _ '^0 


1 ° = 


'3 m° = 


'2 rh° = ^ 


V2 


( 22 ) 


V 2 ’ ^ V 2 ’ V2 

the related complex null tetrad. A basis of self-dual two forms is [5] 

Uab = 2'fh[alb], Vab=2k[amb], lTo& = 2 (rnfaWf,] -f ) (23) 

(note that the only non-zero contractions are UabV°^ = 2 and WabW°^ = —4.) These can be used to expand 


1(7 = + 4'i(C/W -f WU) -f 4'2(DC/ + UV A WW) + 4'3(DW -f WV) + ^liVV 

where UV stands for UabVcd, etc., and 

^0 := Cabcdk°m'>k°m\ := Cabcdk°l^k°m<^ 

^2 := Cabcdk°m^m°l\ ^3 := Cabcdk°l'^m°f 
^4 := Cabcdm°l^fh°l‘^. 

is the complex span of {e“, e^, eg} and, in this basis, 

-f 4-2 - 54'4 5(«'4-4 'o) 

Q"j=\ 1(^-4 - 4-0) + d'2 + 5^'4 l(«'l + ^-s) 

y - ^3 i(^3 + ^3) -2^12 J 


(24) 


(25) 


(26) 
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If Q has three distinct eigenvalues, then the algebraic type of the spacetime is I. If instead two of them equal, say 
Ai = A2 =: A, the space is of type II if dim(ker(Q — AI)) = 1 or type D if dim(ker(Q — AI)) = 2. Finally, in the 
case in which all of three eigenvalues are identical (then necessarily equal to zero, since = 0), the Petrov type will 
be be III, N or O, if dim(ker(Q)) = 1, 2 or 3, respectively. The matrix (26) representing Q'^j can be put in normal 
form (that is, diagonal or Jordan form) by acting on it with elements of SO{3,C). As SO{3,C) is isomorphic to 
SO{3, 1)^ (c.f. equation (16)), any such transformation is uniquely associated with a (proper, orthochronous) Lorentz 
transformation on the spacetime, and this transformation, acting on the null tetrad above, produces what is called a 
principal tetrad. The transformation leading to the normal form of Q, and therefore the principal tetrad, is uniquely 
determined in the cases of Petrov type I, II and III (neither nor P of the unique principal tetrad give a PND in 
type I spacetimes [6]). For a type D space, however, there is a 2-dimensional residual 17(1) x IR>o^ = subgroup 
of S'0(3,1)7 ~ SO{?>, C) of boost and rotations preserving the normal form (and thus the PNDs): 

—>■ ak°', 1“ —>■ m“ —>■ m“ —>■ (27) 

In this case, and and 1“ are aligned along the repetad PNDs, i.e., they satisfy (6). Principal null tetrad components 
of tensors are said to carry spin-weight s and boost-weight q if under (27) they pick up a factor (e.g-, ^'3 has 

s = q = —1). Truly scalar fields, such as Q+ oc 5i(vl'2^), of course, carry zero weights. 


B. Principal null directions 

An alternative approach to the Petrov classification consists of studying the PNDs of the Weyl tensor, i.e., solving 
equation (5), which is equivalent to 


^[e^a]bc[d^f]^ ^ 0 . 

Starting from a generic null tetrad with associated Weyl scalars (25) we find that (see (24)) 


(28) 


2^[e^a]bc[d^f]^ ^ (^ 9 ) 

so the k°‘ vector of the tetrad is a PND if an only the 4/0 component of the Weyl tensor in this tetrad vanishes If we 
apply a null rotation (boost) to the given null tetrad around 

la ya ^ 

k^ -k zzr + zm^ + zm‘^, 

^ -k 2^, (30) 

the resulting will sweep the set of null directions as z moves in the complex plane, avoiding only the 1“ direction, 
wich corresponds to z = oo {S'^ = complex plane plus point at infinity). So we can calculate 5 'q(z) in the primed 
tetrad and solve the fourth-order equation 5 'q(z) = 0 to find the four PNDs. It can easily be checked using (25) that 

^'f){z) = ^'o + 4z5'i -k 6 z ^^'2 + 4z^^' 3 -k (31) 

Generically (type I spaces), there will be four different solutions Zj,j = I, 2, 3, 4 corresponding to four PNDs. The 
special cases are those for which the polynomial (31) has repeated roots, and can be classified according to the 
partitions of 4 as 


Petrov type 

PNDs 

I 

{1111} 

II 

{211} 

D 

{22} 

III 

{31} 

N 

{4} 

0 

{-} 


TABLE I: Petrov type according to the multiplicity of the roots of the polynomial in (31). 

In particular, for type D there are two double roots, and the corresponding will satisfy the stronger equation (6). 
Type O corresponds to conformally flat spaces Cabcd = 0. 
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III. LINEAR PERTURBATIONS 


Let (/ab(e) be a monoparametric family of vacuum solutions with gab(O) =: gab of type D. Assume e“(e) is an 
orthonormal tetrad of the metric gabi^), smooth in e, and such that the associated null tetrad ( 22 ) has A:“( 0 ) and 
?“(0) aligned along the two repeated PNDs of the type D background gab, be., they satisfy ( 6 ). If A(e) is a curve in 
SO{3, l)t with A(0) the identity, then the tetrad e^(e) := A“/ 3 (e)e“(e) satisfy this same condition (this is sometimes 
called the “tetrad-gauge ambiguity”). In any case 4'o(e) = (I/4)C'(e)U(e)U(e) (using (24) and an obvious notation) 
and 


<5«'o := 


_d 

de 


4'o = jiSC VV + C6VV + CVSV) 

e=0 ^ 


\6C VV, 


(32) 


since C5VV = CVSV = 'I' 2 i 5 I 4 f,I 4 “^ = 0. Equation (32) implies that iJ'I'o is tetrad-gauge invariant. The reader can 
check that i 5 'I '2 and are also tetrad-gauge invariant, with 


= \5C UU. (33) 

Note from (23), (27), (32) and (33) that iJ'I'o (( 54 ' 4 ) has spin weight two and boost weight two (minus two and minus 
two respectively). 


A. Perturbed eigenvalues 


To hrst order in e: 


= 'S2 + e6'i>2, 4'i(e) = e,54'„ i = 0,l,3,4 (34) 

Inserting this in (26) we find that the perturbed eigenvalues to order e are 

Ai(e) = 4-2 -f (( 55'2 - \/( 54 'o(I^' 4 )e, (35) 

A 2 (e) = ^'2 + (^4'2 + \/(5To(5^'4)e, (36) 

A 3 (e) = -24'2-2d4'2e, (37) 

and the eigenvalue splitting ( 8 ) is 

2e a/(5To(5^'4, (38) 


the branch choice of the (complex) square root being responsible of the sign ambiguity anticipated in ( 8 ). It is 
important to emphasize that (JTq and d 4'4 are booth free of the tetrad-gauge ambiguity and that they carry opposite 
spin and boost weights (see the discussion around equation (27)). Thus (38) is a well defined scalar field that carries 
information on the distortion of the Weyl tensor due to the perturbation, information that is missing, e.g., by the 
perturbed curvature scalars in the even sector of the Schwarzschild perturbations. 

If either d'l'o = 0 or d 4 t 4 = 0, the space degenerates into a type D or II, depending on the dimension of the eigenspace 
ker(Q — A 2 I). We will comment on these algebraically special perturbations in section IV C. 


B. Splitting of the PNDs 


Replacing (34) in (31) gives 

P{z) := -f 4d4'iez -f 6(^'2 -f ed 4 ' 2 )z^ -f AS'^s.ez^ + SAtiCz'^. (39) 

The equation to be solved is P{z) = 0 up to order e. Note, however, that the solutions z± of the simpler equation 
0 = d'I'oe + 64'2Z^, 








(40) 


satisfy 


P(z±) = 0 + O(e3/2), 


(41) 
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i.e., they are (to order e) two of the four solutions of P{z) = 0. Since —>■ 0 with e, these two solutions are easily 
guessed to be those related to the splitting of A:“ into two different PNDs. Explicitly, in the dominant order we have 


kUe) 






(42) 


According to the discussion between equations (22) and (31) the other two solutions of P{z) = 0 should be near the 
unperturbed repeated PND P, which corresponds to z = oo in = CU{oo}, thus we expect the other two solutions 
to behave as an inverse power or e (c.f. [9]). To obtain these, we can either switch to a; = 1/z, or work with null 
rotations around P and solve the equation '£'4 = 0. In either case we arrive at 


llie) := P±e^/^ 



(43) 


It is important to note that (42) has zero spin weight and boost weight one, and thus define a PND, for which the 
overall scaling is irrelevant. Similarly (43) carries zero spin weight and boost weight minus one. The non-analytical 
character of the splitting, discussed in some detail in [9], can be avoided by a re-parametrization of gab{e). 


IV. APPLICATIONS 

A. Gravitational perturbations of the Schwarzschild black hole 

For the Schwarzschild solution 


dr 


ds^ = —f dt^ + {dO"^ + sin^ 0 dtp'^) , / = 1 — 


2 M 


we use the orthonormal tetrad 


=5 = A( 


61 = 


6n — 


r sin( 


A. 

d(j> 


(44) 


(45) 


then A:“ and P in (22) are along the repeated PNDs. The perturbed metric admits a series expansion using a basis 
constructed from harmonic tensors for [2, 3], labeled (£, m, ±), which can be obtained by applying a differential 
operator to the standard spherical harmonic scalars and are classihed into even (+, or scalar) and odd (-, or 

vector) types according to their behavior under the discrete parity isometry (0,(p) —>■ {tt — 9, (p + tt) The linearized 
Einstein’s equations reduce to two-dimensional wave equations for functions called the Regge- 

Wheeler and Zerilli equations. In terms of these potentials we have found that 


<54/0= ^ A±^(t,r)y2^'’™)(0,<).), 
54'4= ^ B±^(t,r)ri^2 


-2 




(46) 

(47) 


where. 


A 


(£,m) 


3iM l{e + 2)\ 

"S^V (£- 2 )! 


- 1 


(M — r) 


d(j)~ 

dr 


r — 2 M / dt dr"^ 


'(f+ 2 )! 

(£- 2 )! 


2M — r d'^cj)^ 

2 r 2 


1 d'^4) 


Kir) d(j) 


L{r 


, d'^cj) 
dtdr 

dd^ 




2r drdt 2r^D(r) dr 


B 


3iM (£ + 2)1 


{i,m) 


8r3 




{£- 2)1 

2M-r d‘^(j)+ 


.d<f) /r — 3M\ d(j) 

+’'U-2M)ai +’'(2"-’-)^ 


2r2(r — 2M)D{r) dt 
2 d‘^4> 


'{£ + 2)\ 

{£-2)\ [ 2 r 2 


dr'z 




1 d'^<j)'^ K{r) d<jP 


L{r) 


dtdr 

dd'^ 


-f 


N{r) 

4r^D{r) 

f£i£+l) 3M 

2 ~ 

N{r) ^ 


(48) 


2r drdt 2r^D{r) dr 2r‘^{r — 2M)D{r) dt Ar'^D{r) 


(49) 



stands for g^j.g |^]^g normalized spin weight s spherical harmonics on [2], and 

^(r) = (£ + 2)(£-l)r + 6M (50) 

Js:(r) = (^ + 2)(£-l)(M-r)r-6M2, (51) 

L{r) = {£ + 2){e-l){SM-r)r + 6M^, (52) 

N{r) = {£ + 2f(l + !)£{£ - 1)^^ + 6M{i + 2f{l - ifr"^ + + 2){l - l)r + 12M^. (53) 


Note that the eigenvalue splitting (38), as well as the repeated PNDs splittings (42) and (43), being proportional to 
and/or V'f'4 will contain multiple harmonics even if the metric perturbation contains a single non-zero 


B. Gravitational perturbations of the Kerr black hole 

Teukolsky equations [10] are a set of separable partial differential equations for linear fields on type-D backgrounds. 
Two of them describe the behavior of il'I'o and for the type-D background (e.g., a perturbed Kerr black hole), 
assuming a background null tetrad with and 1°' aligned along the repeated PNDs. Although the connection of 
theses quantities with the corresponding metric perturbation is rather intricate [11], Teukolsky equations -unlike the 
Regge-Wheeler and Zerilli equations for a Schwarzschild background- are particularly well suited to our purposes since 
they give the quantities needed in (38), (43) and (42). 

It was shown in [12] that for well behaved (meaning satisfying suitable boundary conditions at the horizon and infinity) 
non-stationary black hole perturbations and (5'I'4 uniquely determine each other. In particular, = 0 if and 
only if 54*4 = 0, and this corresponds to a trivial perturbation. In view of (42) and (43), both repeated PNDs split 
and the perturbed metric is type I. Non-stationary perturbations are those relevant to the black hole stability issue 
and are the ones we focus on in this work. 


C. Chandrasekhar algebraic special modes and Schwarzschild’s naked singularity instability 


In his 1984 paper [13], Chandrasekhar dealt with the problem of finding perturbations of black holes in the Kerr- 
Newman family such that one of (l'I'o,54'4 vanishes while the other does not. It follows from the comments in the 
previous subsection that these perturbations do not satisfy appropriate boundary conditions at the horizon or at 
infinity of a Kerr-Newman black hole. The requirement that d'l'o = 0 or 54*4 = 0 leads to an algebraic condition 
(the vanishing of the Starobinsky constant) that gives a relation among the black hole parameters, the harmonic 
number of the perturbation and its frequency uj (perturbations behave as ~ for pure modes). Although Chan¬ 
drasekhar’s algebraic special (AS) modes in the Kerr-Newman family diverge either at infinity or at the horizon of 
a black hole background, for naked singularities in the Kerr-Newman family some AS are indeed very relevant, as 
they satisfy appropriate boundary conditions both at infinity and at the singularity, and they grow exponentially 
with time (i.e., have a purely imaginary w). The existence of these modes was indeed crucial to prove the linear 
instability of the negative mass Schwarzschild solution and of the the super-extreme Reissner-Nordstrdm solution (for 
the super-extreme Kerr solution, however, none of the AS modes satisfy appropriate boundary conditions and other 
methods were required to establish its linear instability). 


The instability of the Schwarzschild solution (44) with M < 0 is due to the existence of a family of even/scalar 
solutions of the Zerilli equation of the form 

_i_ r(r — 2M)^ 




{t + 2){£-l)r + £>M 


exp 


k{r — t) 
2M 


k = 


(£ + 2 )! 

6 (^- 2 )!' 


(54) 


These were found in [14], then recognized in [15] to agree the the AS modes in [13]. The facts that they behave 
properly for r G (0,oo) (keep in mind that M < 0 and that £ >2 for non-stationary perturbations, for further details 
see [14]) and grow exponentially with time, signals the instability of the naked singularity. For the perturbations (54) 
we find 


< 54-0 = 0 , 


Jc 

<5'I^4 = ^(r-2M) 


k-1 


i£ + 2){£-l) 


6M 


exp 


( k{r-t) \ 
\ 2M ) 


^-2 


(55) 


According to (42) and (43) the double PND will remain double whereas /“ will split, the perturbed spacetime 
being being of type II according to Table I. Note, however, that this can only be accomplished by fine tuning the 
perturbation to restrict to the modes (54), a generic perturbation will also contain the stable, oscillating modes, and 
the PNDs will split into four, that is, to type I. 
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V. DISCUSSION 

We have found explicit formulas for the splitting of the repeated PNDs of type D spacetimes under gravitational 
perturbations, and also for the splitting of the repeated eigenvalue in (4). These are given in equations (42) and (43) 
and (38). These are observable (gauge invariant) effects of the perturbation on the background geometry that do not 
require higher than two derivatives of the metric, in contrast to differential invariants. 

In view of (42), (43) and the results in [12], perturbed black holes within the Kerr-Newman family suffer a PND 
splitting to type I, except for stationary perturbations, which by the black hole uniqueness theorems are restricted to 
changes in the mass and/or angular momentum parameters, and therefore trivially keep the type D structure. This 
gives sense to the notion that the “tangent vector” hab of a curve gab{f) at a black hole metric go&(0) will “stick 
out” of the set of algebraically special metrics. Note, however, that boundary conditions play a crucial role in these 
assertions: the example in Section IV C shows the flow of the type D Schwarzschild naked singularity to a type II 
spacetime, triggered by the instability. This flow, however, can only be occur for fine tuned initial conditions allowing 
only the (infinitely many) modes (54). A generic perturbation will contain modes other than these and will therefore 
split the two repeated PNDs into four single PNDs. 

As a final comment, the non-analytical behavior of the PNDs in the perturbation parameter (the dominant order in 
the perturbed PNDs is e^/^) is to be expected from the polynomial character of the PND equation and the confluence 
of the solutions as e —>■ 0. This fact was clarified in [9], whose results are in total agreement with ours. 
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